Irregular Dynamics in a One-Dimensional Bose System 
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We study many-body quantum dynamics of J-interacting bosons confined in a one-dimensional 
ring. Main attention is payed to the transition from the mean-field to Tonks-Girardeau regime using 
an approach developed in the theory of interacting particles. We analyze, both analytically and 
numerically, how the Shannon entropy of the wavefunction and the momentum distribution depend 
on time for a weak and strong interactions. We show that the transition from regular (quasi-periodic) 
to irregular ("chaotic") dynamics coincides with the onset of the Tonks-Girardeau regime. In the 
latter regime the momentum distribution of the system reveals a statistical relaxation to a steady 
state distribution. The transition can be observed experimentally by studying the interference 
fringes obtained after releasing the trap and letting the boson system expand ballistically. 

PACS numbers: 05.45Pq, 05.45Mt, 03.67,Lx 



The study of the quantum many-body dynamics of 
a Bose-Einstein condensate (BEC) is a major challenge 
in the physics of trapped gases. The mean- field Gross- 
Pitaevskii framework has successfully described an im- 
pressive set of experimental data pj, but it fails when 
genuine quantum many-body (QMB) correlations are im- 
portant (e.g., in quantum phase transitions QUI, in the 
design of quantum information devices Q and interfer- 
ometers jl| , and near dynamically unstable regions where 
quantum corrections appear on a logarithmically short 
time-scale 

QMB effects also become particularly striking in low- 
dimensional systems, where the competition between 
quantum fluctuations and statistical properties is par- 
ticularly enhanced. Recent experimental achievements 
in effective one-dimensional (ID) harmonically confined 
quantum degenerate systems have stimulated several 
efforts to understand their main properties §]. The ID 
regime can be achieved in optical/magnetic traps when 
the radial degrees of freedom are frozen by tight trans- 
verse confinement. The Hamiltonian of the boson gas is 
given by the Lieb-Liniger model Q, where the two-body 
interaction is assumed to be point-like. The thermody- 
namical properties of this Hamiltonian, and its excita- 
tion spectra have been calculated analytically by Licb 
and Liniger [9(- An unexpected feature, predicted by Gi- 
rardeau [ifj, is the onset of fermionization when n/g — > 
(here n is the particle density, and g is the interatomic 
coupling constant which is inversely proportional to the 
ID interatomic scattering length |ll|)- In this Tonks- 
Girardeau (TG) regime, the density of the interact- 
ing bosons becomes identical to that of non-interacting 
fermions (while, of course, the wave-function keeps the 
bosonic symmetry). On the other hand, in the opposite 
limit, n/g — > oo, the system is described in the mean- 
field (MF) approximation as a weakly interacting boson 



gas. The crossover between these two regimes occurs near 
n/g ~ 1. In this region the MF approach breaks down 
and more complicated two-body correlations become cru- 
cial. 

While the thermodynamics of a ID boson gas is fairly 
well understood, systematic studies of many-body dy- 
namical properties of this system have been started only 
recently (see |l'2l and references therein). In this con- 
text, one should mention investigation of general quan- 
tum (beyond MF) properties of the dynamics of trapped 
bosons, which can be considered as the next frontier in 
BEC studies. In pursuit of this goal, in this Letter we 
study numerically the quantum dynamics of N interact- 
ing Bose particles on a one-dimensional ring of length L. 
Our main result is the discovery of a dynamical transi- 
tion from regular to (quantum) irregular dynamics when 
n/g — > 1. This transition is accompained by a linear in- 
crease of the Shannon entropy of the wave packets, that 
is directly related to the vanishing of the interference 
fringes which occurs after releasing the confinement and 
letting the bosons expand ballistically. 

Our model is specified by the Hamiltonian, 
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Here — ata^ is the occupation number operator; 
a\ and a s are the creation-annihilation operators; and 
efc = Air 2 k 2 /L 2 . The classical analog of this model, which 
is originated from the nonlinear Shrodinger equation has 
been analyzed in other papers. In the context of our 
study, the most interesting results have been obtained in 
|l2j | where the interchange between regular and irregu- 
lar dynamics has been explored (without the connection 
with the Tonks-Girardeau regime). Note, also, that our 
model Q is purely quantum and the Hamiltonian is writ- 
ten in the operator form (see also ^3). Specifically, the 



occupation numbers for single-particle levels are not as- 
sumed apriori to be very large. To analyze this model, 
below we use the methods developed for quantum sys- 
tems of interacting particles with complex behavior (see, 
e.g. 0] and references therein). 

In what follows, we explore the situation when all 
bosons initially occupy the single-particle level with the 
angular momentum k = or, equivalently, at t = the 
system is in the unperturbed (g — 0) ground state \^f g ). 
Note that the total momentum is conserved in time. Our 
main interest is the evolution of the system for different 
values of the control parameter n/g. This type of exper- 
iment can be realized since the ratio n/g can be manip- 
ulated by tuning the interatomic scattering length. 

It is convenient to label single-particle states \k) ac- 
cording to their momentum k — 0,±1,±2... ± m, .... 
Then, any many-body state \j) can be represented as 
I-.., n J _ m , n J Q , n J m , ...), where n J k represents the num- 
ber of particles in the single-particle level characterized 
by the momentum k of the j-th state. In our numerical 
studies we need to take a finite number of particles TV 
in a finite number £ — 2m + 1 of single-particle states. 
Clearly, the number N of atoms and i of levels should be 
chosen in a consistent way in order to have results which 
can be extrapolated to an arbitrary number of atoms. In 
a ID geometry on a ring, N particles define the small- 
est spacing, which corresponds to the largest value of the 
momentum m f» N. In our numerical simulations, we 
satisfy the latter relation when choosing changing values 
of N and m. 

The transition from the MF to the TG regime can be 
understood using the following handwaving argument. 
Having all particles initially in the lowest state with 
k = 0, let us estimate the strength of the interaction 
necessary to move two particles from the unperturbed 
ground state to the upper (and lower) single particle 
levels k = ±m. One can expect that this interaction 
will result in an ergodic filling (in time) of all single- 
particle states. The energy required for the shift is ap- 
proximately m 2 / L 2 w N 2 /L 2 , and the matrix element 
of the interaction between the corresponding states is 
< j\V\f >~ g^fN(N - Vj/L. Equating these values 
one can obtain n — N/L ~ g that is associated with the 
crossover from the bosonic to fermionic regime. Dynam- 
ically this crossover is reflected by a rapid depletion of 
the fraction of particles in the lowest single-particle state 
with k = 0. Note that at n/g ~ 1, the ratio Nq/N ~ h, 
with N = (do* do). 

One of the most appropriate quantities to characterize 
the dynamical properties of this system is the Shannon 
entropy of the wave packet, 

s(t)=-x;i* i (t)i a ]iii* i (t)i a (2) 

3 

Here &j(t) — (j\$f(t)) is the projection of the wave func- 
tion onto the noninteracting many-body basis. After 
switching on the interaction, the wave function evolves 
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FIG. 1: Entropy as a function of the rescaled time gt for fixed 
m — 6 and different values of n/g. Data are given for N = 
6. The saturation of S(t) for small values of n/g manifests 
the onset of the Tonks-Girardeau regime. The dashed line 
corresponds to the theoretical prediction J3J. 

according to the total Hamiltonian, ^f(t) = e~ lHt '$> g (0) 
and spreads over the unperturbed basis. 

The Shannon entropy, defined by Eq.© allows one 
to estimate the effective number of unperturbed many- 
body states that are involved in the dynamics due to 
the inter-particle interaction, N e ff(t) w exp{S'(t)}. A 
few examples of the time dependence of the entropy are 
reported in FigQ] One can see that for n/g ^ 1 the 
entropy oscillates in a regular way, while for n/g <C I 
there is a generic linear time dependence for a short time 
followed by a saturation at later times. 

While saturation is due to the finiteness of our Hilbert 
space, the linear growth in time of the entropy was shown 
[l4| to be associated with the onset of chaos and thermal- 
ization in close systems of randomly interacting particles 
(see details in [l5| and 0). Specifically, for a strong 
enough interaction between particles the linear increase 
of S(t) is due to an exponential increase in time of the 
number of many-body basis states involved in the dy- 
namics. Although our model is not random, off-diagonal 
matrix elements determined by the particle interaction 
strongly fluctuate in time due to the different number of 
particles occupying single particle levels. This fact allows 
us to apply the results obtained in [Tij for our model. 

For the case when the probability of remaining in the 
originally excited state decreases exponentially, Wo(t) = 
exp(— Ft), it was found that the entropy increases lin- 
early, 

S(t)=TthiM, (3) 

apart from a generic quadratic increase that occurs at 
small times. Here the value of T is determined by the de- 
cay of the probability to stay in the unperturbed ground 
state, and M = m is the number of many-body unper- 
turbed states directly coupled to the initial state. As was 
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FIG. 2: Rescaled entropy as a function of n/g for N = 6 
particles and different m as indicated in the caption. 




FIG. 3: Particle density in the momentum representation for 
four values of the interaction n/g, and for different rescaled 
times gt = 0, 0.2, 0.4, 0.6, 0.8 (lines are the thicker for later 
times). 



shown in j 1 31 ] , the linear growth of entropy appears in the 
regime for which T is proportional to perturbation g. 

We have carefully checked the predictions for the time- 
dependence of S(t) given in Ref.0]. I n particular, we 
have confirmed that the number of excited basis states 
grows in time exponentially fast before saturation. The 
saturation occurs because the total number of many- 
body states that we had in our numerical simulation is 
finite. We also have found that the probability Wo (t) de- 
cays exponentially in time. As a result, we were able to 
verify the increase of entropy given by Eq.Q and indi- 
cated as a dashed line in Fig^ 

We summarize our main results in Fig|21 Here S aV e is 
the mean value of the entropy averaged over time after 
the saturation, and Njj is the dimension of the Hilbcrt 
space (the total number of many-body states defined by 
m and N). Entropy rescaling has been done in order 
to avoid its trivial dependence on the size of the Hilbert 
space. In this case the plateau for n/g <C 1 is indepen- 
dent of the choice of the number of particles and the 
number of levels. One can see that the renormalized en- 
tropy is a quantity that is very sensitive to the transition 
MF-TG regime. Also, as one can see, slight variations in 
to as N value does not change the results. This indicates 
the universal character of the entropy in the transition 
region. Note that the limiting value of the renormalized 
entropy is not 1, that would correspond to the prediction 
of random matrix theory. This deviation clearly indicates 
the influence of the dynamical nature of the Hamiltonian. 

We now discuss how our predictions can be tested 
experimentally. Since the entropy of wave functions 
is not an easy quantity to measure experimentally, we 
need to relate it to an observable which can be mea- 
sured directly. In cylindrical coordinates the boson field 
operator (defined on a ring of radius R in the plane 
2 = 0) has the form, ip(p, 9, z) = 8{z)5{p — R)ip(9). Here 
the angle-dependent part can be expressed in terms of 



plane waves, ip(9) = X^^fe^ • Fourier transforming, 
x(p) — J d?,r {jj(p, 9, z) e~ l P' r , and taking into account 
that p — pe y and p ■ r = ppsm9, we have, 

X(p) = Cj2 & k[ dOe lkm e- ipRs ' m{9) (4) 

where C is a normalization constant. Using periodic 
boundary conditions, k = ^-n = with 77 = 0, ±1±2...., 
one obtains 

+00 

X(P) = C JnipR)a n (5) 

n— — 00 

where J n (x) is the Bessel function of order 77. 

Expresssing the wave function at time t in terms of 
the basis states \j) = \n J _ o:i , ...n J , ...n 3 +00 ), we obtain 
an expression for the occupation number distribution 
n(p,t) = {^(t)]^ (p)x(p)\^p(t)) . This quantity gives the 
probability density for finding a boson with momentum 
p along the y axis at time t. Using our expression we 
obtain, 

+00 

j,f Z,s=-oo 

(6) 

Since the total momentum is zero, the selection rules, 
(J\a,a s \j') = njSjj'61 s hold, where 77j is the number of 
particles occupying the single-particle level with momen- 
tum / in the j-th basis state. Finally, we obtain, 

+00 

n(p,t) = CY\Ut)\ 2 E J l(P R ) n i CO 

j i = -oo 

where C _1 = iV _1 J dp n(p,t) . 
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FIG. 4: Visibility as a function of time, for four values of 
rescaled interaction n/g. Here is N = 6 and m = 6. 

Results are reported in FigO When n/g >> 1, the 
momentum distribution does not practically change in 
time. All particles remain mainly in the condensate dur- 
ing the evolution. Close to the transition point n/g < 1, 
the distribution begins to be flat, without relaxing to 
an asymptotic distribution. In this intermediate region, 
the occupation number distributions reveal a mixture of 
very slow relaxation and strong oscillations in time that 
was first observed |13| in a similar model and termed as 
"non-ergodic behaviour of interacting bosons." 

Contrary to this, for n/g -C 1, one can see a clear 
relaxation of the occupation number distribution to a 
steady state distribution. In fact, the onset of relaxation 
to the steady distribution means that the behavior of the 
model can be well described by statistical methods 0] ■ 

In order to characterize qualitatively the onset of re- 
laxation in the TG regime, we propose that one observe 
experimentally the fringe visibility, defined as 

^(^) \Imax I min\ / ij- max ~i~ ^rain) (8) 

Here we chose I ma x = n(p = 0, t), l m i n = n(pg, t) where 



Po ~ 2.41 is the first zero of Jo- See FigQ] This quantity 
allows one to distinguish between periodic and aperiodic 
behavior of the occupation number distribution. These 
different time behaviors are shown in Fig0| 

As one can see, the TG regime is characterized by a 
fast decay of the visibility with subsequent fluctuations 
that persist in time. These fluctuations are additional 
evidence of the statistical nature of the dynamics when 
the interparticle interaction is strong. When n/g 3> 1 
(MF regime), the visibility oscillates, revealing collapses 
and revivals of the condensate fraction population, as the 
entropy does in the same regime. See FigQ] Therefore, 
there is a strict correspondence between the dynamical 
evolution of the entropy and the visibility of the momen- 
tum distribution. The intcrfcrnce properties of a thermal 
Tonks gas localized on one side of a ring has been studied 
in 0. 

In conclusion, we have studied the dynamics of a Bose- 
Einstein condensate model on a torus by paying atten- 
tion to the role of the interparticle interaction. Specif- 
ically, we considered the situation in which initially the 
ground state only, defined in the absence of the interac- 
tion, is populated. By switching on the interaction we 
have found different regimes depending on the strength 
of interaction. The first regime (MF regime, n/g 3> 1) 
is characterized by regular dynamics associated with pe- 
riodic oscillations in time of the Shannon entropy Eq. 
and of the fringe visibility Eq. (JSJ of particle density. In 
the Tonks- Girardeau regime, n/g <C 1, the Shannon en- 
tropy grows linearly in time in accordance with our ana- 
lytical estimates. Even if these estimates are obtained for 
completely random models, they give remarkable agree- 
ment with our data. Moreover, this approach allows one 
to understand the mechanism of the transition from the 
mean field to Tonks gas regime. 
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